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Abstract. 

It is known that two-dimensional lower branch Tollmien-Schlichting waves described 
by triple-deck theory are always stable for planar supersonic flows. Here the possible 
occurrence of axisymmetric unstable modes in the supersonic flow around an axisymmetric 
body is investigated. In particular flows around bodies with typical radii comparable with 
the thickness of the upper deck are considered. It is shown that such unstable modes exist 
below a critical nondimensional radius of the body ao. At values of the radius above ao 
all the modes are stable whilst if unstable modes exist they are found to occur in pairs. 
The interaction of these modes in the nonlinear regime is investigated using a weakly 
nonlinear approach and it is found that, dependent on the frequencies of the imposed 
Tollmien-Schlichting waves, either of the modes can be set up. 


This work was supported under the National Aeronautics and Space Administration 
under NASA Contract No.NASl-18107 while the authors were in residence at ICASE,NASA 
Langley Research Center, Hampton VA23665. 


1. Introduction. 


It is now well-known that lower branch Tollmien-Schlichting waves in incompressible 
boundary layers are governed by triple-deck theory. The linear description of Tollmien- 
Schlichting waves near the lower branch of the neutral curve in the wavenumber-Reynolds 
number plane (in the limit of large Reynolds number) was given by Smith (l) whilst the 
extension into the corresponding weakly nonlinear regime was given by Smith (2). It was 
shown in the latter paper that nonlinear effects stabilize the linear growth of the waves so 
that as they propagate downstream a local equilibrium state is reached. Both two- and 
three- dimensional modes can be described by this approach and it was shown by Hall and 
Smith (3) that when nonlinear effects first become important nonparallel effects cannot be 
ignored. Further downstream the waves adjust in a quasi-parallel manner to the local flow 
conditions. 

As the disturbance moves downstream its amplitude grows and ultimately the large 
amplitude high frequency state of Smith and Burggraf (4) is set up. Here the disturbance 
has a Stokes layer at the wall and resonant triad interactions can occur (Smith and Stewart 

( 5 )). 

For supersonic planar flows the complex eigenrelation for two-dimensional Tollmien- 
Schlichting waves has only stable solutions. Here the possible existence of unstable ax- 
isymmetric modes in supersonic flows is discussed. We show that such modes exist when 
the nondimensional radius of the body about which the fluid flows falls below a critical 
size. This is consistent with the planar result which is achieved in the limit of large body 
radius. 

In fact we find that unstable modes occur in pairs so that in the finite amplitude 
regime a nonlinear interaction problem must be set up in order to see which is the pre- 
ferred mode. This is done using the approach of Hall and Smith (3) and we find that 
the interaction is characterized by a pair of coupled cubic nonlinear ordinary differential 
amplitude equations. The solutions of these equations depend on the values of the coeffi- 
cients appearing in them. The latter coefficients depend on the radius of the body so that 
a different bifurcation structure is found as the radius varies. We show that, depending on 
the radius and frequencies of the Tollmien-Schlichting waves, two types of stable solution 
are possible. 

In the next section we formulate the triple-deck problem which governs lower branch 
Tollmien-Schlichting waves in supersonic flows. In §3 these equations are solved for dis- 
turbances so small that linear theory is a valid approximation. In §4 the weakly nonlinear 
nonparallel interaction of the two possible instability modes of §3 is discussed. In §5 we 
discuss the results of the interaction problem of §4 and draw some conclusions. 

2. Formulation 

In this study we are concerned with the stability of an axisymmetric boundary layer 
which forms on a cylindrical body of radius a*, in a uniform supersonic stream of velocity 
Uoo directed along the axis of the cylinder. 
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Suppose that L* denotes a typical streamwise lengthscale, v*oo the kinematic viscosity 
of the fluid , then the Reynolds number Re is defined to be 


Re = UooL'/v^, 


( 2 . 1 ) 


this will be taken to be large throughout this paper, and so the parameter 


e = Re~ l/% 


( 2 . 2 ) 


may be taken to be small. 

It will be assumed that 

a = ^ = 0(D (2-3) 

which defines the scale of the radius of the body (this follows the order of a * considered 
by Kluwick et al (6). 

The scalings of the various components of the solution then fall broadly into the 
triple-deck model. The following non-dimensional variables are defined 


X = 


P = 


v = 


X 


£ 3 L* 

P* ~ P*oo 

PloUg ’ 

u* ’ 

^oo 


X — 1 

I 5 ”’ 

U = 


r = 


u 

U*~' 

^oo 


C = 


U2o' 


f U'oot* 
e*L* 


e*L*' 


(2.4) 


Here x * , r* denote streamwise and radial co-ordinates respectively, u*,v* the corresponding 
velocity components, c* the speed of sound, the fluid density in the external flow, 
P*(Poo) the pressure (in the external flow), and t* denotes time. 

The stability of the flow is investigated at a downstream location. Here the boundary 
layer thickness is 0(e 4 L*), which is somewhat thinner than the radius of the body; we are 
thus concerned with streamwise locations relatively close to the leading edge of the body. 
At the axial position in question, we suppose that the skin friction is A e~ A /L*, and in 
Section 4 we shall incorporate non-parallel effects into our study by allowing streamwise 
variations of A. 

The first layer of the triple-deck to be considered is the upper deck, wherein F = 0(1). 
The solution develops (Kluwick et al (6)) as follows 


u = 1 + e 2 ui(X,f,t)+, 
v = <r 2 vi(X,f,t)+, 
p = e 2 p 1 (X,r,t)+, 
c = Af" x + e 2 ci(X,r, t)+, 


(2.5) 
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where Moo is the Mach number of the external flow. The problem may be reduced to 
the solution for the pressure, which is governed by the axisymmetric form of the Prandtl 
Glauert equation, viz 

(1 - Ml^Ptxx + rPir + Piff ~ °* ( 2 * 6 ) 

The solution to this will be discussed later. 

The main deck corresponds to the transverse scale 


r — a 


= 0 ( 1 ), 


and is the same as planar flow to leading order, with the solution taking the form 

u = Uo[y) + eA[X,t)U{ ) (y)+, 
v = - e 2 Ax{X , t)U 0 (y)+, 

P = Ro(y) + eA{X,i)Ro y (y), 

p = e*P{X,t) + - 

Uo(y) and Ro(y) are the axial velocity and density distributions respectively, in the undis- 
turbed boundary layer. 

Taking the limit of the v expansion as y — ► oo in (2.8) yields an inner boundary 
condition for (2.6) namely 

Pif 1^=3— (2.9) 

whilst if y — ► 0, the no-slip condition on the wall is violated, and this demands the presence 
of a lower deck wherein 

y = y/e = 0(l), 

and 

u = \Y +eU{X,Y,i)+, 

v = e 2 V{X,Y,t)+, 

_ (210 > 
p = .Ro(O) 4 - spi , (X,Y , t) 4 - . 

At this stage it is convenient to scale out a number of the physical parameters. Following 
Kluwick et al (6) , this is achieved with the following transformations 

X = CI|M^-1|-I(T W /T 00 )IX, 


X = Cb 

Y = C* 
P = Ci 

U = C* 

V = c§ 
A = C* 
a= C* 


Ml -l 
Ml -l 
Ml -l 
Ml- 1 
Ml -l 


Ml- 1 


{T w /Too)*Y, 

-ip, 

{TJToo^U, 
* (Tui/Too)^, 
{T w /Too)*A, 
(TjT^a, 
(T w /Too)t. 


( 2 . 11 ) 
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Here C is the Chapman constant which arises from the linear viscosity law 


(**Ko = C(T/T„), 


( 2 . 12 ) 


and the nondimensional wall temperature is denoted by T w . The governing equations in 
the lower deck reduce to the planar incompressible form at leading order, namely 


dU_ dU dU dHJ_ 

dt +U dX + V dY ~ dY 2 
dU dV n 

dX dY 
dP n 

dY ~ ’ 


dP_ 

dX' 


(2.13) 


with boundary conditions 

U = V - 0 on y = 0, 

U -► X[Y + A{X, t)] as Y -»• oo. (2.14) 

The problem is closed by means of solving (2.6) subject to (2.9), giving (symbolically) 

P = £(A). (2.15) 


In general, the system is nonlinear, and requires a numerical treatment; however, 
progress may be made by seeking solutions which are slight perturbations about the undis- 
turbed state U = Ay, V = P = A = 0, and this will form the basis for the following section. 

3. Linear Stability 

Here the basic flow is slightly perturbed by a small amount h « 1, and the solution 

is thus written 

U = \Y + hU 1 + ---, 

V = hV x + • • • , 

P = hP 1 + ---, 

A — hA\ + • • • . 

A solution is sought of the form 



U\ — U\{Y)E\ + c.c., 
V 1 =V l (Y)E l +e.c. i 
Pi = P\E\ + c.c., 

Ai = A\E\ + c.c.. 


(3.2) 


Here 


E\ = exp[t(aX — fit)], 


(3.3) 
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and ’c.c.’ denotes a complex conjugate, Ui{Y) and Vi(y) are complex functions of Y 
alone, and Pi,Ai,a and fi are, in general, complex constants. 

The problem for Ui{Y) then reduces to 

i(aY - U)Uiy = Uiyyy , 

with 

Ui = 0 on y = 0 

A A 

Ui — *• A\ as y — > oo. 

At this stage it is necessary to determine the solution of (2.6), subject to (2.9) which 
yields: 

^ _ iaKo{icta)Ai 
Pl ~ ~K\ (taa) - 

This solution (see Ward (7) and Kluwick et al (6)) is chosen such that information prop- 
agates downstream along the characteristics of (2.6) . K n (z) denotes the modified Bessel 
function of order n argument z. 

The system comprising (3.4), (3.5) is closed by utilizing (3.6) with 



(3.4) 

(3.5) 


VlYY Iy=0= 1 , 


(3.7) 


which arises from evaluating (2.13) on Y = 0. 

A 

The solution for I7iy may be written down in terms of Airy functions, and a non- trivial 
solution is only possible if the following dispersion relationship is satisfied 


-A 3 At'o = (mA)*jr 0 ( taa) 


where 

= -(*a)^n/a Xt 
Ai' 0 = At'(eo). 


(3.8) 


(3.9) 


It is well known that neutrally stable Tollmien-Schlichting waves arise from (3.8) when 
Imag{a} = Imag{ fi} = 0. It is also worth noting at this stage that in the limit of a — ♦ oo, 
i.e. the flat plate limit, (a fixed), (3.6) degenerates to the plane supersonic interaction 
condition 

Pi tail, (3.10) 

which is known to exhibit only decaying solutions (see for example Ryzhov and Zhuk (8) 
). On the other hand, in the limit of a very thin cylinder (or needle) , namely o->0, then 

Pi ~ — a 2 aii log a, 


(3.11) 
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which is akin to the interaction condition found in certain jet and channel flows, (see 
Smith ( 9 )), which does admit growing Tollmien-Schlichting instabilities (Bogdanova and 
Ryzhov ( 10 )). From these somewhat qualitative observations, we may suspect that prior 
to a thorough investigation of the roots of (3.8) , growing Tollmien-Schlichting instabilities 
are possible (at least over a range of a). 

A search was carried out (using Newton interation) for neutral roots of (3.8), i.e. for 
which Imag{a} = Imag{fl} = 0 , with A = 1 (indeed, A may here be scaled from the 
problem, but is retained for the purposes of the following section). No roots of this type 
were found for a > ao = 0.003 (approximately, although it is likely that an infinite number 
of exponentially decaying roots do exist in this regime). For a < ao it was found that two 
neutrally stable roots exist for each value of a. Figure 1 shows the variation of neutral 
frequency with a, whilst Figure 2 shows the variation of wavenumber with a. 


Notice that both fl and a become unbounded at a — ► 0 , and this can be confirmed 
analytically. First consider the lower branch, marked L in Figures 1 and 2 . We seek an 
asymptotic solution to ( 3 . 8 ) in the limit as a — ► 0 with £o = 0 ( 1 ) , and a = 0 (a -1 ). Then 
it follows that 

Ko[iaa) . . in v . „ . 

y ~ -taa(ln(a) + — + 7 + . (3.12) 

where -7 = 0.577215 . . ., and so 


% ~ j + -]• 


(3.13) 


Hence the leading order term on the right-hand-side of this expression has argument 7 r/ 6 , 
and it is well known from planar incompressible stability studies (see Smith ( 1 ), Hall and 
Smith (3)) that the right-hand-side of (3.13) has this same argument when 


£0 - -2.298*^, 


which corresponds to 


A»'(6>) 

Mm 


. 1 

— t ^ . 


At leading order (3.13) then becomes 

aia | log(aa) |~ A^. 


(3.14) 


(3.15) 


(3.16) 


This is a transcendental equation linking a to a (see Duck ( 11 ) for a similar type of 
situation); however it may be solved by successive approximations for a » 1 . 

The first approximation gives 

a = 0(a^). 


(3.17) 
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The second approximation gives 


a 


=( 


7\i 


4 a I log a 




(3.18) 


(this process may be continued further, with the inclusion of terms like log | log a | etc., 
but we shall terminate the approximation at this stage). 

Taking (3.18), we obtain an asymptotic expression to fi, namely 


n 


2.298 ( 


7Al 


4a I log a 


f) T - 


(3.19) 


The asymptotic results (3.18) and (3.19) are indicated on Figures 1 and 2 by broken 
lines, and the agreement with the computed results as a — ► 0 is seen to be satisfactory. 
(We cannot expect better agreement in view of the likely largeness of the correction terms 
to (3.18), (3.19) which have been omitted.) 

We now turn to consider the asymptotic form of the upper branch (denoted by U on 

Figures 1 and 2) as a — ► 0. Since we are concerned with neutral modes, it is convenient to 

set 

Co = ~i*R, (3.20) 

R = n/oci, (3.21) 

and to seek solutions to (3.8) for which R =| Co |— 1 ► oo, a — *■ 0, a = 0(a -1 ). With these 

restrictions, (3.8) becomes 


i.e. 


—A* (ta)^[ln(aa) + ^ + If *i ] 

~ \ 3 [-i*R + e~ 2 ^R~^] i 


(3.22) 


a^a[/n(aa) + y + 7H ]~A£[-iE + e *SF R ]. (3.23) 

We now take the leading order real and imaginary parts of this equation separately to 
yield two (real) equations linking a, i2(ft) and a, namely 

a$a | ln(aa) |~ \$R, (3.24) 

«*•! ~ R~*‘ (3.25) 

The leading order terms in the solution for this system are 


a~2 ^A?7^a * n ^(|ina|) ^ [l + ^ na — ] 

Vl " L 28 | Ina | J 

fl ~ Ina |)»[l - 1 1 ^ 11 ]. 


(3.26) 
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These asymptotic results are shown as ’dot-dashed’ lines on Figures 1 and 2. The 
agreement is not as satisfactory as that for the lower branch asymptotic behaviour consid- 
ered previously; here the presence of the many log, log log etc. terms considerably degrades 
the convergence of the asymptotic series. However results obtained by the authors (not 
shown here) for a « 10 -4 , indicate the proposed upper branch structure to be correct. 
The existence of two modes of instability, found above, suggests that interaction between 
them must be taken into account when studying the nonlinear problem. We consider this 
aspect in the following section. 

4. Wave interactions in the weakly nonlinear regime 

Here we shall discuss the nonlinear nonparallel interaction of two axisymmetric 
Tollmien-Schlichting waves at a value of a less than the cut-off value Oo. The interaction is 
essentially similar to that given by Hall and Smith (3) for three-dimensional interactions 
in incompressible flow so all of the details will not be repeated here. We first note that the 
Tollmien-Schlichting waves vary on the triple-deck length scale e~ 3 in the linear regime. 
Now we are concerned with their possible equilibration on a longer length scale centred 
on the downstream position where they are locally neutrally stable. It was shown by Hall 
and Smith that the appropriate length scale where nonparallel effects are most important 
is 0(e*). At distances large compared to e* downstream of the neutral position the 
disturbances suffer no nonparallel effects to the order given in our calculation. 


The linearized form of the evolution equation for a Tollmien-Schlichting wave of am- 
plitude B in an O(frl) neighbourhood of the neutral position x n can be derived from (3.8) 
by writing 

V • • , $ & 

X = a — , ta = ta + — — , 

el dX 


and expanding to give 


ajj ~ 

+ eXB = 0, 
dX 


(4.1) 


where 


1 aK i aKp _ 2t £ Atp . 2 i tl At 0 » 

3 ?a Ko (taa) Ki (taa) 3 a ^ 0 ' 3 VQ ctAt^ 


(4.2) 


where 


Ai 0 = Ai(£ 0 ),Xo = f Ai(()d( 
U 0 


Thus on the basis of linear theory it follows that in an 0(e!) neighbourhood of x n , B is 
given by 

_ . , 

B = constant e - 3 (4.3) 


This linear growth will take the size of the initial disturbance to a size where nonlinearity 
cannot be neglected. This occurs when the x velocity component in the lower deck is of 
order 
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Suppose then that at z = z R the two modes with (a, fi) = (ay, fly), j = 1,2 are 
neutrally stable. We define Ej , j = 1,2 by 


_ . r a,(x — x n ) ^ , 

Ej = exp t{ - - ~o Hyt}, 


(4.4) 


where fcy and fly are now real quantities. We must expand ay and fty in the form 


ay — ayo + £c*yi + • • • , 
fly = flyo + sflyi + • • • , 


(4.5) 


where (ayo, flyo) and (ayj,flyi) are chosen such that the waves are neutrally stable at 
x n . Thus (ayo, fly,,) are as given in §3 whilst (fcyi, flyi) can be found by continuing the 
expansion procedure of that section to next order. For the purposes of this calculation we 
do not need to know (ayi, Qyi) so we give no further details of their determination here. 
We suppose that the waves have grown sufficiently for nonlinear effects to be important 
so that in the lower deck the appropriate expansion is 

(u,v,p) =(eAY,0,0) + {4[e*Cfii,e^Vi 1 ,e^P 11 ][l + 0(e)]£ 1 

+ BleZUi 2f e?V 12 ,e*P 12 ]ll + 0(e)]E 2 + A 2 letu 21 ,e*V 21 ,e$P 21 ]ll + 0(e)] 

+ B 2 [e*U 22 , e$V 22 , eip 22 \\l + 0(e)]£?f + AB[e*U 24 ,e$V 24 ,e 7 *P 24 ][l +0(e)] 

+ AB[e? U 22 , V 23 , e^p 23 ][i + 0(e)]£i^2} + c *c. 

+ I A | 2 [e%U 2 z, 0, 0] [l + 0(e)]+ | B | 2 [e^l726)0>0] + ••• 

(4.6) 

Here the factor [1+ 0(e)] is needed because ay and fly are ’neutral’ to order e and the terms 
represented by • • • include the effect of the interactions between first harmonic, mean and 
fundamental terms. These interactions produce terms proportional to E\ and E 2 and it is 
the solvability condition on the equations for the coefficients of these quantities which give 
the required amplitude equations for A(X) and B(X). We further expand A in the form 

A = A(x n ) + et^lx + • • • , (4.7) 

to account for the variation in the shear stress away from the neutral position . In the 
main and upper decks similar expansions can be written down but there is no interaction 
between the modes so that the solutions of the equations are essentially unchanged from the 
linear analysis of §3. For that reason we shall not write down the expansions in these layers 
but merely quote the required matching conditions for the lower deck where needed. It is 
a routine matter to substitute the expansions (4.6), (4.7) into the Navier-Stokes equations 
and equate like powers of e for each Fourier component. 

The zeroth order equations for the fundamentals show that (C/ 11 , Vn,Pu), 
(Ui 2 ,V 12 ,P 12 ) are simply the eigenfunctions of the linear problem appropriate to (a, fl) = 
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(c*io,ftio) and (a, fl) = (0:20,^20) respectively. At next order the equations for the func- 
tions proportional to Ei,E 2 i EiE 2 t EiE 2 are obtained. We find that (t^2i>^2i>-P2i) sat- 
isfies 


{— 2tflio + 2taioA(x n )y}£/2i + TjjiA = — 2ia.ioP2\ + U 2i — {*«io^ii + ^ii^n}> 
2 tQ!l 0^21 + ^21 = 

^21=721 = 0^ = 0, 

( 4 . 8 o, 6 , c) 

and the required matching conditions imposed by the main deck are 

U21 -* A 21 X(x n ),Y - oo; 2 ai 0 A 2 i = • ( 4 - 9 ) 

The function {U 22 ,V 22 , P22), {U 2Z , V 23 ,P 23 ) and {U 24 ,V 24 ,P 24 ) satisfy similar inhomoge- 
neous differential systems with the wavenumber and frequency changed appropriately and 
the nonlinear term on the right hand side of ( 4 . 8 a) replaced by —{ia. 2 oU \ 2 + V\ 2 U[ 2 }, 
— {*(0:10 + oc 2 o)UnVi 2 + Vi 2 U' ll + Vu U' 12 } and — {*(<*10 — o. 2 o)UixU\ 2 -\-Vi 2 U [ 1 + VuU' 12 } 
respectively. In the absence of the inhomogeneous term the above systems would of course 
not have a solution. Finally at this order the mean flow correction functions U 2 s and U 2 q 
are found to satisfy 

K » = 7 n^i + VixV[ 1, II" = V 12 U [ 2 + V 12 U[ 2 , 

U 2 & = U 2 q = 0 , y = 0 ; U 25 ,U 26 —+ 0 , Y — ► 00. 


The next order system in the expansion procedure leads to equations which determine 
an, <221 the correction terms in the expansions of the neutral wavenumbers. However the 
system obtained next involves terms proportional to E\ and E 2 produced by nonlinear 
interactions and nonparallel effects. Thus, unless the required orthogonality condition is 
satisfied, there will not be a solution of the equations obtained at this order. After some 
manipulation the required solvability conditions yield 


dA 

dX 

dB 

dX 


= -ci XA + axA | A | 2 +h A | B | 2 , 
= —c 2 XB + a 2 B\B | 2 +b 2 B | A | 2 . 


( 4 . 10 a, 6 ) 
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(4.11o-/) 


Here c\ and c 2 are given by ( 4 . 2 ) with ( 0 ,a) = (U 10,0:10) and (n,a) = (n20>«2o) respec- 
tively whilst 01,61,02,62 are defined by 

<■1 = J°° T l (Y)F[(Y)dY/d l , 6, = r r i (Y)Fi(Y)ir/iu 

d t = { 0 -r l( 0)^2)[2. ai o + + $!2^)]£M<*>) 

iiri(taioa) v K 0 (*aio<») Ki[iai 0 a) n 

+ n{Y\(x n )U n )'ri{Y)dY }, 

Jo 

a 2 = r r 2 {Y)G\{Y) dY/d 2 , b 2 = [ °° r 2 {Y)G’ 2 {Y) dY/d 2 , 

Jo Jo 

. ( m\ K o(*a20«) ro . . 2 (K\{ia 2 oa) K[{ict 20 a h,_ r . . 

d 2 = {-»- 2 ( 0 )Tr-r £ 2 *«20 + 0 = 20 a (irh £ + ITT- £) ^12(00) 

Ki\ia 2 oa) s Ko{ta 20 a) K 1 (ta 20 a) ,i v ’ 

+ [°° {Y\(x n )U l 2 )'r 2 (Y) dY. 

Jo 

Here Fi , F 2 , G\ , G 2 are given by 

■f’l = — {*o:io(^24^ii + U 2 iUn) + Fnl^ + V 2 \U[ X + V n U 21 }, 

F 2 = —{iot\oU\ 2 U 2 4 + i&\oU 2 QU\\ + %OL\oU 22 ,tj\ 2 + V\ 2 U 2 ^ 

+ V 23 u [ 2 + V 2 A U [ 2 + V u UU + V X 2 U' 2Z }, 

Gi = —{ia 2 o(U 2 QUi 2 + U 22 U \ 2 ) + Vi 2 U 2Q + V 22 U [ 2 + Vi 2 U 22 }, 

G 2 = —{iot 2 oUiiU 2 4 + * 0120 ^ 25^12 + * 0 ^ 20 ^ 23^11 + V11O24 
+ V 23 ^1 + ? 2A U’ n + V\\U 2Z + V l2 U’„}, 

whilst ri(y) and r 2 (Y) are the adjoint eigenfunctions defined by 


(4.12a, 6, c,d) 


0 0 


d ( Pi 

~£Y ( $3 | = — ( ^ ® t'flyo + \(x n )iajoY 


r i 


0 1 


0 

.2 



, , a vn -f^o ( to: ,n a) r ,( 0) 

qj = 0, on y = 0; ry = 0, asF -> 00 , py(oo) = -J - — r — -, j = 1,2. (4.13) 

AJ\i\%OLjoa) 

5 . Results and Discussion 

Following Hall and Smith (3) it is an easy matter to generalize (4.10) to the case 
where the wave has moved a distance 0(h 2 ) where h is small compared to unity but large 
compared to e*. Here the variable X is now defined by 


- _ (x - x n ) 


X = 


h\ 
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7 

and the size of the disturbances discussed in §4 are now changed from, for example, u ~ e* 
in the lower deck to u~ eh. We then find that (4.10) becomes 


dA 

dX 

dB 

dX 


— — Ci AA + o,\A | A | 2 ~\-biA | B | 2 , 
= -c 2 AB + a 2 B\B | 2 +b 2 B | A | 2 , 


(5.1a, b) 


where A = A(x n ) + Ah 2 H . A further generalization of (5.1) is obtained by detuning the 

frequency of either of the Tollmien- Schlichting waves by an amount sufficient to introduce 
further linear terms in these equations. Without any loss of generality we suppose that 
the ’A’ mode is determined so that (5.1) becomes 


dA 

dX 

dB 

dX 


— —c\AA + gi A + a\A \ A | 2 -\-b\A \ B | 2 , 
= -c 2 AB + a 2 B | B | 2 + b 2 B \ A | 2 . 


(5.2) 


Here g\ is a complex constant with g\ r > 0 if the A mode is the most unstable and g\ r < 0 
otherwise. In order to discuss the solution of (5.1) we write 


«=MI 2 . 1>=\B\\ 


so that 

1 dS . \ . . 

= (— ClrA + gir)6 + + girty, 

\% , ^ 

r— = -c 2r Arj) + a 2r il> + b 2r 6tf>. 

* dX 

Apart from the trivial solution 6 = t(> = 0 the following finite amplitude solutions are 
possible. 

a. 6 — 0 xp = c 2r A(L 2 r 

b. 6 = (cir A - d Xr )aXr , 0 = 0, 

C. 6 = [a 2 r(ci r A — dir) — C 2r bi r ]\ai r a 2 r — 6lr62r] _1 , 

0 = [C2r - Sb 2r }a 2 r. 

Here a subscript r denotes the real part of a quantity. We refer to a, b, c as the ’pure 
B\ ’pure A’ and the mixed mode respectively. The range of values of A for which the 
different solutions exist and their stability properties depend on the coefficients in (5.3). 
This problem has been discussed by Keener (12) for the case when ci r , c 2r , ai r , a 2r , b\ r and 
b 2r are all negative. 

The results in Table I show that f 21 - < ^ ^ for all the values of a used in the 

Dir Ci r air 

calculations. It is straightforward to show that in this case whichever mode bifurcates first 
remains stable at all values of A and suffers no secondary bifurcation to the mixed mode. 
The second ’pure’ mode to bifurcate is initially unstable but then undergoes secondary 
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bifurcation to the mixed mode. The mixed mode is always unstable but the pure mode 
from which it bifurcates is stable after the bifurcation. Thus at sufficiently high values of A 
either pure mode is a possible stable state. However if the modes are set up by increasing 
A we would of course expect to observe just the most unstable linear mode. These results 
are illustrated in Figure 3 for the case where the pure A mode bifurcates first. When a 
tends to its critical value from above the two modes become identical so that for example 
a\ — * 62 and 2a\ — ► 61 . This served as a useful check on the calculations. 
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Figure 1. The variation of the neutral wavenumber with 
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Figure 2. The variation of the neutral frequency with a. 
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TABLE I 


a 

a /b 
r r 

b /a 
r r 

c /c 
r 

.0028 

.590 

1 .684 

.935 

.0026 

.592 

1.672 

.916 

.0024 

.592 

1.666 

.907 

.0022 

.589 

1.660 

.903 

.0020 

.585 

1.655 

.903 

.0018 

.580 

1.651 

.905 

.0016 

.574 

1.647 

.904 

.0014 

.567 

1.645 

.915 

.0012 

.559 

1.643 

.922 

.0010 

.549 

1.642 

.931 

.0008 

.537 

1.643 

.943 

.0006 

.526 

1 .642 

.952 
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